Low temperature transport in granular metals 
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We investigate transport in a granular metallic system at large tunneling conductance between 
the grains, gr 3> 1- We show that at low temperatures, T < grS, where 5 is the single mean energy 
level spacing in a grain, the coherent electron motion at large distances dominates the physics, 
contrary to the high temperature (T > grS) behavior where conductivity is controlled by the 
scales of the order of the grain size. The conductivity of one and two dimensional granular metals, 
in the low temperature regime, decays with decreasing temperature in the same manner as that 
in homogeneous disordered metals, indicating thus an insulating behavior. However, even in this 
temperature regime the granular structure remains important and there is an additional contribution 
to conductivity coming from short distances. Due to this contribution the metal-insulator transition 
in three dimensions occurs at the value of tunnel conductance g^ = (1/6tt) \n(Ec/5), where Ec is 
the charging energy of an isolated grain, and not at the generally expected g% oc 1. Corrections 
to the density of states of granular metals due to the electron-electron interaction are calculated. 
Our results compare favorably with the logarithmic dependence of resistivity in the high-T c cuprate 
superconductors indicating that these materials may have a granular structure. 

PACS numbers: 73.23Hk, 73.22Lp, 71.30.+h 



A great deal of research in the current mesoscopic 
physics focuses on understanding properties of granular 
metals (see P,0>I3)- The interest is motivated by the 
fact that while their properties are generic for a wealth of 
strongly correlated systems with disorder, granular met- 
als offer a unique experimentally accessible tunable sys- 
tem where both the interaction strength and degree of 
disorder can be controlled. 

The key phenomenon revealing the most of the under- 
lying physics is transport, where the effects of interac- 
tions play a crucial role. The processes of electron tun- 
neling from grain to grain that govern electron transfer, 
are accompanied by charging the grains involved after 
each electron hop to another grain. This may lead to a 
Coulomb blockade, and one justly expects this effect to be 
of the prime importance at least in the limit of weak cou- 
pling. It makes it thus clear, on a qualitative level, that 
it is the interplay between the the grain-to-grain coupling 
and the electron-electron Coulomb interaction that con- 
trols transport properties of granular metals; yet, despite 
the significant efforts expended, a quantitative theory of 
transport in metallic granular systems is still lacking. 

A step towards formulation such a theory was made 
recently in Ref. 0. It was shown that depending on the 
dimensionless tunneling conductance gx one observes ei- 
ther exponential-, at <C 1, or logarithmic, at gx 3> 1 
temperature dependence of conductivity. The consider- 
ation in Ref. y was based on the approach developed 
by Ambegaokar, Eckern and Schon (AES) 5] for tunnel 
junctions. This technique however, as shown in Ref. 
applies only at temperatures T > grS, where S is the 
mean energy level spacing in a single grain, in this regime 
the electron coherence does not extend beyond the grain 



size. At low temperature region, T < gxS, the effects of 
the electron coherent motion at distances much exceeding 
the single grain size a must be included, thus this impor- 
tant regime is not described by the AES approach . 

Although experimentally the low temperature regime 
is well within the experimental reach [jj, |2| , it has never 
been addressed theoretically so far. The important ques- 
tion whether the system is a metal or becomes an in- 
sulator, in other words, whether the conductivity of the 
granular metals at large conductances remain finite in 
the limit of T — > is still open. 

In this Letter we investigate the low-temperature con- 
ductivity of granular samples focusing on the case of large 
tunneling conductance between the grains, > 1. To 
this end we develop a technique that goes beyond the 
AES approach and includes effects of coherent electron 
motion at distances larger than the size of the grain. 
Without the Coulomb interaction the granular system 
would be a good metal in the limit, gx 3> 1, and our 
task is to include the charging effects into the theory. 
We find that at temperatures, T < gxS properties of the 
granular metal depend on the dimensionality of the array, 
and corrections to the conductivity and density of states 
due to the effects of Coulomb interaction are similar to 
those obtained in Ref. for a homogeneous metal. Thus 
at low temperatures the systems behaves essentially as a 
homogeneous metal contrasting the case of large temper- 
atures, T ^> grS considered in Ref. 0. 

This in particular means that at large conductances the 
3D system is a good metal. On the other hand, at <?t •C 1 
a granular sample is in the insulating state. Therefore 
a 3D system should exhibit a metal-insulator transition 
at the critical value of the conductance gx, such that 
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samples with conductances qt > g^ are metals and their 
conductivity remains finite at T — > while samples with 
9t < 9t are insulators and their conductivity vanishes 
at T -> 0. 

The main results of our work are as follows: (i) We 
find the critical value g% of the tunnel conductance at 
which the metal-insulator transition in 3D occurs 
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(1/6tt) HEc/8), 



(1) 



where Ec is the charging energy of an isolated grain, (ii) 
We find the expression for the conductivity of a granular 
metal that includes corrections due to Coulomb inter- 
action and holds for all temperatures as long as these 
corrections are small. The corresponding answer can be 
conveniently written separating the correction due to the 
contribution from the large energy scales e > gxS from 
that coming from the low energy scales e < grS. Denot- 
ing corrections as 8<ji and 8<ii respectively we have 



a = <tq + Sai + 5<J2, 



(2a) 



where Co = 2e 2 f/Ta 2_d , with a being the size of the single 
grain is the classical Drude conductivity for a granular 
metal (spin included). Correction Sai in Eq. (J23) con- 
tains the dimensionality of the array d only as a coeffi- 
cient and is given by the following expression 0], 



1 



2ndgj 



In 



9tE c 



max (T, g T S) 



(2b) 



On the contrary the correction 5<J2 in Eq- l|2al) that is im- 
portant only at temperatures T < Sgr strongly depends 
on the dimensionality of the array 




D = 3, 
D = 2, 
D = l. 



(2c) 



Here a = / °° dx 
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Dth(x) + x/sinh 2 (x)] « 1.83 



and /3 = J °° dxx~ 3 ^ 2 [coth(x) — xj sinh 2 (izi)] w 3.13 are 
the numerical constants. For a 3D granular system a 
temperature independent term of the order 1 / gj> has been 
subtracted in the first line in Eq. Q2c[). 

Corrections Sai and 5cr 2 are of a different origin: the 
correction Sai comes from the large energy scales, e > 
gxS where the granular structure of the array dominates 
the physics. The fact that this correction is essentially 
independent of the dimentionality d means that the tun- 
neling of electrons with energies e > grS can be consid- 
ered as incoherent. On the other hand, correction 5a2 
in Eq. (|2"c|) is similar to that obtained for homogeneous 
metals long ago jfj and comes from the low energy scales, 
£ < grS, where the coherent electron motion on the scales 
larger than the grain size a dominates the physics. 

It is important to note that in the low temperature 
regime all temperate dependence of conductivity comes 



from the correction (5cr 2 . At the same time, in this regime 
the correction Sai, though being temperature indepen- 
dent, still exists and can be even larger than 8a2- 

When deriving Eqs. (J2J) we neglected possible weak lo- 
calization corrections that may originate from quantum 
interference of electron waves. This approximation is le- 
gitimate if a magnetic field is applied as in Ref. or 
dephasing is strong due to inelastic processes. 

Now we turn to the description of our model and the 
derivation of Eqs. (0): We consider a d— dimensional ar- 
ray of metallic grains with the Coulomb interaction be- 
tween electrons. The motion of electrons inside the grains 
is diffusive and they can tunnel from grain to grain. In 
principle, the grains can be clean such that electrons scat- 
ter mainly on grain surfaces. We assume that the sample 
in the absence of the Coulomb interaction would be a 
good metal. For large tunneling conductance we may also 
neglect the nonperturbative charging effects (discretness 
of the electron charge), which give an exponentially small 
(as exp(— #(?t)) contribution to the conductivity. Al- 
though we assume that the dimcnsionless tunneling con- 
ductance <7t is large, it should be still smaller than the 
grain conductance, go , such that gr < go- This inequal- 
ity means that the granular structure is still important 
and the main contribution to the macroscopic resistivity 
comes from the contacts between the grains. 

The system of weakly coupled metallic grains can be 
described by the Hamiltonian 

^[^ t (r l )^(r,)+^ t (r 3 )^(r 4 )], (3) 



H — Ho + H c 
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where is the tunneling matrix element corresponding 
to the points of contact and rj of i-th and j— th grains. 
The Hamiltonian Hq in Eq. describes noninteracting 
isolated disordered grains. The term H c describes the 
Coulomb interaction inside and between the grains. It 
has the following form 



(4) 



where C'ij is the capacitance matrix and is the operator 
of electrons number in the z-th grain. In the regime un- 
der consideration one can neglect the coordinate depen- 
dence of a single grain diffusion propagator. The electron 
hopping between the grains can be included using the di- 
agrammatic technique developed in Refs. 0, 0, which we 
outline below. 

Electron motion in a random impurity potential within 
a single grain can be considered using the standard dia- 
grammatic techniques described, for example in Ref. [3]. 
Electron hoppings between the grains can be considered 
in a similar way assuming that tunneling matrix ele- 
ments between the grains are random variables obeying 
the Gaussian statistics and correlated as 

{tki,k 3 tk 3 ,ki) = ^ (^k lt k s 4 2 ,fe 4 + Ski,k 4 $k 2 ,k 3 )' (^) 



3 



-®- 



<8^ 



a) 



b) 



(a) 



_L 



FIG. 1: Self energy of the election Green function averaged 
over impurity potential inside the grains and over tunneling 
elements between the grains. Averaging over the impurity 
potential is represented by the dotted line (a) while tunneling 
elements are represented by crossed circles (b). 



where t is related with the avarage intergranular conduc- 
tance as gx = 2nt 2 /5 2 . The average Green function is 
defined by the Dyson equation where self energy, shown 
on Fig. 1 has two contributions: The first contribution 
(a) corresponds to scattering inside a single grain while 
the second (b) is due to processes of scattering between 
the neighboring grains. Both this processes result in a 
similar contribution ~ sign (a;) to the electron self-energy 
thus on the level of single particle electron Green func- 
tion intergranular scattering results only in small renor- 
malization of the relaxation time r 



2dg T S, 



(6) 



where To is the electron mean free time in a single grain. 

The next step is to consider the diffusion motion of 
electron through a granular metal: Diffusion motion in- 
side a single grain is given by the usual ladder diagram 
that results in the diffusion propagator 
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where coordinate dependence was neglected since we as- 
sume the zero dimensional limit for a single grain. Tun- 
neling between the grains is accounted for in a similar 
way, such that the total diffusion propagator is given by 
the ladder diagrams shown on Fig 2a. This results in the 
following expression: 
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T |fi|+fe q ' 
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where e q = 2gy^ a (l — cosqa) with a being the lat- 
tice vectors. For small quasimomenta q <C a -1 we have 
E q — > grSa 2 q 2 such that the propagator (JSJ describes the 
diffusion motion on the scales much larger than a with 
effective diffusion coeeficient D = gra 2 5. 

The same ladder diagrams describe the dressing of in- 
teraction vertex as it shown on Fig. 2b. The dressed 
vertex can be used to obtain the polarization operator, 
that defines effective dynamically screened Coulomb in- 
teraction (Fig 2c): 



V(n,q) = 



C(q) 
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FIG. 2: These diagrams represent (a) Dyson equation for dif- 
fusion propagator (b) interaction vertex dressed by impurity 
and intergranular scattering (c) Screened Coulomb interac- 
tion. 



The conductivity of the granular metals is given by the 
analytical continuation of the Matsubara current-current 
correlator. In the absence of the electron-electron inter- 
action the conductivity is represented by the diagram (a) 
in Fig. 3 that results in high temperature (Drude) con- 
ductivity do which is defined below Eq. I|2al) . First order 
interaction corrections to the conductivity are given by 
the diagrams (b-e) in Fig. 3. These diagrams are analo- 
gous to ones considered in Ref. @ for the correction to 
the conductivity of homogeneous metals. We consider the 
contributions from diagrams (b,c) and (d,e) separately: 
The sum of the diagrams (b,c) results in the following 
correction to the conductivity 

cr 27 



-J— ImV / 

Indgx J 



dw7(w)e q y(w,q). (10) 



where 7(w) = ^jwcoth^, and the potential V(ui,q) is 
the ana;itic continuation of the Screened Coulomb poten- 
tial with dressed interaction vertices included attached at 
both ends 

2£c(q) 



(eqS - faS) (4e q £'c(q) - iu) ' 



(11) 



The above expression was simplifed using that the charg- 
ing energy Ec(q) = e 2 /2C(q), expressed in terms of 
the Fourier transform of the capacitance matrix C(q) 
is much larger than 6. Performing the integration over 
the frequency and summing over the quasimomentum q 
in Eq. I|1U|I with the logarithmic accuracy we obtain the 
correction (|2b|l . One can see from Eq. I|10[l that the con- 
tribution 8<j\ in Eq. I|2bjl comes from the large energy 
scales, e > gr5 such that at low temperatures the loga- 
rithm is cut off on the energy scale grS. 

To obtain the total correction to the conductivity of 
granular metal the two other diagrams, (d) and (e) in 
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FIG. 3: Diagrams describing the conductivity of granular 
metals: the diagram (a) corresponds to <7q in Eq. 12aH and 
it is the analog of Drude conductivity. Diagrams (b)-(e) de- 
scribing first order correction to the conductivity of granular 
metals due to electron-electron interaction. The solid lines 
denote the propagator of electrons and the dashed lines de- 
scribe effective screened electron-electron propagator. The 
sum of the diagrams (b) and (c) results in the conductivity 
correction Sai in Eq. [|2al . The other two diagrams, (d) and 
(e) result in the correction 5o2- 



Fig. 3 should be taken into account. These diagrams 
result in the following contribution to the conductivity 



5a 2 
cto 



2g T S 



/ dui j(oj) Im 



^(w,q) Ea sin2 (q a ) 



(12) 

In contrast to the contribution 5a\ in Eq. I|ll)[l . the main 
contribution to the sum over the quasimomcntum q in 
Eq. 112fl comes from the low momenta, g < 1/a. In this 
regime the capacitance matrix, C(q) in Eqs. and 
(|12|1 has the following asymptotic form 
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Using Eqs. (|1 1113(1 . we obtain the result for the correction 
Sa 2 in Eq. (|2c[l . This correction has a physical meaning 
similar to that of the Altshuler-Aronov correction 6] de- 
rived for homogeneous disordered metals. 

Comparing our results in Eqs. f2f with those obtained 
in Ref.U using the AES functional we see that the correc- 
tion to the conductivity obtained in Ref. is equivalent 
to the correction 5o\ in Eq. (|2al) . which corresponds in 
our approach to the sum of diagrams (b) and (c) in Fig. 3. 
The correction Sa 2 in Eq. (12al) becomes important only at 
low temperatures, T < gxS where AES functional is not 
applicable. While in our approach both corrections to the 
conductivity must be small 5<Ji, Sa 2 <C (To the method of 
Ref.H gives a possibility to show that for T gxS the 
dependence of the conductivity is logarithmic so long as 
a/e 2 a 2 ~ d > 1. 

It follows from Eq. (|2"c|> that at low temperatures, 
T < gxS, for a 3D granular array, there are no essen- 
tial corrections to the conductivity coming from the low 
energies since the correction 5a 2 is always small. This 
means that the result for the renormalized conductance, 
gx of Ref. (see also ) for 3D samples within the log- 
arithmic accuracy can be written in the following form 



9t{T) = g T - — In 

D7T 



g T E, 



c 



max (gxS, T) 



(14) 



such that it is valid for all temperatures as long as 
the renormalized conductance, gT 3> 1- One can 
see from Eq. l(T4"|) that for bare conductance, gx 3> 
(l/67r) hi(gxEc /S) the renormalized conductance, gx is 
always large and the system remains metallic down 
to zero temperatures. In the opposite limit gx < 
(1 / 6n) ln(grEc / S) , the system flows when decreasing the 
temperature to the strong coupling regime, gx ~ 1 that 
indicates the onset of the insulating phase. We see that 
with the logarithmic accuracy the critical value of the 
conductance gx is given by Eq. (QJ. 

The result for the bare critical conductance in Eq. (JTJ 
agrees with the estimate for g^, that follows from the 
consideration of Coulomb blockade phenomena in a sin- 
gle grain ^l| : the contribution of Coulomb blockade to 
thermodynamic quantities in the regime of strong cou- 
pling is controlled by the factor ~ exjj[—irg(T)], where 
g(T) = gx — (1/Zn) \n(gxEc /T) with Z being the num- 
ber of contacts. Coulomb blockade effects become strong 
at g(T) ~ 1. Taking T ~ gxS and Z — 6 we estimate 
the bare conductance as g^ ~ (l/67r) hx{gxEc /T) that 
coincides with Eq. 

Corrections to the density of states (DOS) can be 
obtained in a simialr way by considering the diagrams 
shown on Fig 4. The diagran (b) results only in the 
energy shift, and therfore is not imprtant, while the dia- 
gram (a) results in the following contribution 
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> Im 

4tt ^ 
q 



tanh[(£-w)/2T] 



(e q (5 - iuj)[eq - wj/4Ec(q)] ' 
(15) 

Here isq is the DOS for noninteracting electrons, e q and 
Ec(<l) were defined below Eqs. IjlUI ) and (|11|) respec- 
tively. Using Eq. (|15fl for a 3D granular array we obtain 
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E c gx 



max (e, gxS) 



(16a) 



where A = g T a 3 J d 3 q /(27r) 3 e~ 1 and e = max{T, e}. 
For £ 3> gxS the correction to the DOS (|16a|l coincides 
with the one obtained in Ref. Q using AES approach. It 
follows from Eq. I|16a|) that for a 3D array of grains, as 
in case with conductivity, the main contribution to the 
DOS comes from the large energy scales, e > gx5. 

Using Eq. i|15|) for a 2D array we obtain the following 
result for the correction to the DOS 



5v 2 



1 



_ e > gz5, 

16g T TT 2 ) In ^ In + 2 In 2 ^ e < g T 6. 



2hx 2 {gxEc/i) 



(16b) 

Using the relation between the tunneling conductance 
and the diffusion coefficient, D = gx& 2 5 one can see 
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FIG. 4: Diagrams representing corrections to the single par- 
ticle density of states 

that the temperature (energy) dependence of the DOS 
for e <§C <7t<5 given by Eq. i|16b(l coincides up to the con- 
stant term with the result for the correction to the DOS 
of the homogeneous metal . 

The logarithmic behavior (|2b|) of the conductivity is in 
a good agreement with experimental findings 0, Q ■ It 
would be very interesting to extend the resistivity mea- 
surements to the low temperature regime where we pre- 
dict the temperature dependence (J2SJ- A similar loga- 
rithmic dependence of resistivity on temperature was re- 
cently found in high-T c compounds La,2-ySr y CuOt± and 
Bj^Sr2- x La x CuOQ+s in a very strong magnetic field 
(lot A possible granularity of these samples was 

suggested in Ref. 0- Recently the microscopic gran- 
ularity was directly experimentally observed in the su- 
perconducting state of Bi 2 Sr2CaCu208+s by the STM 
probe |k| • If we accept that samples studied in 0, 
are indeed microscopically granular, we can compare the 
results of the experiments with our predictions. When 
doing so it is convenient to scale three dimensional con- 
ductivity to the conductivity of CuO planes, a v i ane . Ac- 
cording to our predictions 

da p i ane / dlnT = (e 2 / 1 nfi) fc, (17) 

where the coefficient k — l/2ir in the low temperature- 
and k = 1/d in the high temperature regimes. While in 
the low temperature regime the application of Eq. I|17|) 
is legitimate only under the assumption that electrons 
in different CuO plane are incoherent, in the high tem- 
perature regime the behavior of conductivity according 
to Eq. (|2b|) is logarithmic for any dimension. In this 
regime the real dimensionality d should be replaced by 
d = Z/2, where Z is the (average) number of the contacts 
of each grain with all the adjacent grains. Describing the 
data shown in Fig. 3 of Ref. 0] by our log dependen- 
cies at temperatue T * 5K we extract k ~ 0.4, for Sr 
concentration of y = 0.08 for La,2- y Sr y CuOi [l4j; for 
the Bi2Sr2- x La x CuOe+s compound we find k ~ 0.2 for 
x = 0.84 La concentration, and k ~ 0.3 for x — 0.76. For 
each particular curve the values k extracted from Fig. 3 



of Ref. [12J increase with temperature (especially in case 
of LSCO), this is in a complete agreement with our re- 
sults provided that the "coherent-incoherent" crossover 
occurs at about T ~ 5K . At higher temperatures k no- 
ticeably exceeds 1/2tv, supporting the idea of a granular- 
ity of doped cuprates. 

In conclusion, we have investigated transport proper- 
ties of granular metals at large tunneling conductance 
and obtained corrections to the conductivity, Eqs. (J2J 
and DOS, Eqs. due to electron-electron interaction. 
We have shown that at temperatures, T > grS the gran- 
ular structure of the array dominates the physics. On the 
contrary at temperatures, T < grS the large-scale coher- 
ent electron motion is crucial. Comparison our results 
with experimental data supports the assumption about 
a granular structure of doped high-T c cuprates. 
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